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摘 要 : 在 干旱 少雨 的 山区 开展 小 流域 的 暴雨 山洪 预报 预警 关键 技术 研究 ,对 防 灾 减灾 意义 重大 。 


2018 年 7 月 31 日 新 疆 哈 密 北 部 山区 出 现 特大 洪 惠 


域 大 暴雨 面 雨 量 
资料 (CMPAS) 计 算 了 射 月 沟 水 库 上 游 


量 与 水 利 部 门 事后 调查 数据 较 吻 合 ,最 大 洪 


AT ,发 生 罕见 的 山洪 灾害 ,致使 射 月 沟 流 域 水 库 漫 坝 
溃 口 ,下游 受灾 严重 。 射 月 沟 流 域 气象 观测 站 点 少 且 缺乏 水 文 监测 资料 ,为 客观 定量 分 析 射 月 沟 流 
\ 形 成 的 洪水 汇 水 量 以 及 致 灾 水 库 过 程 。 通 过 采用 空间 插值 法 和 多 源 融 合 逐 时 降水 
雨量 并 进行 检验 分 析 。 根 据 不 同 面 雨 量 驱 动 Floodarea 模 型 
得 出 射 月 沟 水 库 上 游 素 计 汇 水 量 ,结果 表明 :多 源 融合 降水 产品 估算 所 得 最 大 洪峰 流量 和 累计 汇 水 
量 为 1756m.s" ,精确 性 达到 调查 值 的 95%, 射 月 沟 水 


库 上 游 暴 雨山 洪 总 量 为 2.64x107m, 远 超 该 水 库 的 防洪 库容 和 溢洪道 承载 能 力 。 
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件 频 发 ,给 人 民生 命 财产 和 社会 经 济 带 来 极 大 损 
失 ” ,因此 ,认识 地 球 生态 系统 演变 规律 ,并 从 内 
在 机 理 出 发 ,逐步 深入 研究 ,掌握 极端 降雨 引发 的 
暴雨 山洪 变化 规律 ,及 时 规避 洪涝 灾害 和 防 灾 减 灾 
都 非常 重要 。 在 暴雨 山洪 监测 预警 方面 ,西方 发 达 
国家 起 步 较 早 ,欧美 从 20 世 纪 70 年 代 起 陆续 建立 
了 山洪 灾害 预报 预警 系统 ,基本 实现 了 山洪 灾害 的 
实时 监测 和 水 文 模型 的 山洪 预报 预警 5”。 我 国 山 
洪灾 害 防治 起 步 较 晚 , 近 10 a 来 ,我 国 实施 了 山洪 灾 
害 防治 项 目 建设 ,初步 建成 了 山洪 灾害 防御 体系 "， 
并 在 监测 预警 以 及 应 急 避 险 建设 取得 长 足 进 步 ,但 
暴雨 山洪 罕 发 性 强 .成 灾 快 、 破 坏 性 强 .预警 预防 难 
度 大 ,山洪 灾害 仍 是 我 国 洪水 灾害 中 导致 人 员 伤亡 
最 主要 的 灾 种 之 一 。 

目前 ,关于 不 同 尺度 的 洪灾 预警 预报 研究 侧重 
于 峰值 的 预报 、 洪 水 过 程 以 及 面 雨量 的 预报 等 方 
面 。MICHAUD 等 ”通过 中 小 尺度 流域 不 同类 型 水 
文 模型 的 模拟 对 比 表 明 ,降雨 输入 精度 对 模拟 结 
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影响 较 大 ; 张 雪 松 等 "发 现 雨量 站 密度 、 分 布 和 降 
雨 空间 分 布 变 化 均 会 对 水 文 模拟 结果 产生 影响 ; 章 
国 材 等 "基于 GIS 栅 格 数据 的 水 动力 暴雨 洪涝 淹没 
模型 ,对 流域 暴雨 洪水 淹没 过 程 进 行 淹 没 模拟 , 推 
算出 山洪 淹没 不 同等 级 的 风险 雨量 ,取得 了 与 实际 
较为 一 致 的 结果 。 但 是 现 有 的 山洪 灾害 预报 预警 
主要 依据 实时 降雨 和 水 位 信息 ,存在 较 多 空 报 和 预 
和 警 后 来 不 及 避 险 撤离 转移 ” ,尤其 在 山地 小 流域 地 
区 ,在 气象 ,水文 站 点 稀缺 且 分 布 不 均匀 的 情况 下 ， 
开展 短 时 山洪 预警 难度 较 大 ,小 流域 山洪 灾害 动态 
预报 预警 技术 仍然 是 山洪 灾害 防御 急需 解决 的 关 
键 技术 问题 ”“。 因 此 ,通过 现代 气象 监测 手段 和 各 
种 数值 降水 预报 模式 ,在 山洪 灾害 较 严 重 的 地 区 建 
立 小 流域 面 雨量 \ 汇 水 量 等 定量 化 预警 产品 既是 暴 
雨山 洪灾 害 预警 和 灾害 风险 评估 的 技术 关键 所 在 
又 是 防 灾 减 灾 的 迫切 需求 。 

脓 雨 虽然 是 非常 重要 的 灾害 性 天 气 ,但 是 暴雨 
并 不 一 定 会 直接 造成 生命 伤亡 和 财产 损失 ,暴雨 造 
成 的 灾害 是 通过 暴雨 引发 的 山洪 、 泥 石 流 等 次 生 灾 
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害 ERED KRUELA 24 h 累计 降水 达到 50 
mm 时 并 不 会 造成 山洪 灾害 ,但 是 在 北方 地 质 条 件 
和 生态 环境 差 的 地 方 , <50 mm 的 降雨 量 可 能 会 引 
发 地 质 灾害 或 山洪 "。 新 疆 地 形 地 貌 复杂 ,山脉 与 
盆地 相间 排列 ,气候 干旱 , 春 夏季 多 又 十 天气 ,尤其 
阿尔 泰山 北部 .昆仑 山脉 南部 以 及 中 天 山 一 带 的 高 
山 丘 陵 区 多 有 极端 的 暴雨 天 气 发 生 , 易 形成 暴雨 山 
洪 造 成 重大 人 员 伤 亡 和 财产 损失 “'。2018 年 7 月 
29 日 ~8 月 2 日 上 午 (北京 时 间 , 下 同 ) ,哈密 市 北部 
山区 出 现 强 降 水 天 气 过 程 , 尤 其 是 7 月 31 日 01 时 ~ 
14 时 ,哈密 市 伊 州 区 至 伊 吾 县 山区 出 现 大 范围 特大 
俊 雨 ,其 范围 内 射 月 沟 流 域 的 沁 城 乡 小 保 、 射 月 沟 
水 库 及 周边 的 沁 城 乡 3 个 自动 气象 站 分 别 观测 到 
116 mm、105 mm、79 mm 的 强 降水 , 均 突破 有 记录 以 
来 极 值 ;小 保 站 在 31 日 06 时 ~09 时 连续 3h 雨 强 超 
过 全 国 暴雨 小 时 降水 强度 (16 mm ho!) ,降雨 强度 之 
大 在 干旱 少雨 的 西北 地 区 实 为 罕见 (图 1, 图 2), 特 
大 暴雨 造成 小 流域 内 射 月 沟 水 库 漫 坝 溃 口 ,下 游 受 
灾 严 重 。 本 文 就 此 次 过 程 , 研 究 山 区 小 流域 在 气 
象 水文 监测 资料 不 完备 的 情况 下 ,利用 Floodarea 
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图 1 哈密 市 伊 州 区 射 月 沟 流 域 及 周边 气象 站 分 布 图 


Fig. 1 Weather stations of the Sheyuegou Basin, Yizhou 


District, Hami City and surrounding region 


30 上 = 小 堡 站 
= 沁 城 乡 
25 上 = 射 月 沟 水 库 
E 20 上 
ra 15 
Bet 


图 2 MEID AH Yak eA EET N E AI 


Fig. 2 Hourly rainfall of the three automatic rainfall stations 


near Sheyuegou Reservoir on 31st July 2018 


模型 模拟 骏 雨 洪水 致 灾 射 月 沟 水 库 的 水 文 过 程 , 通 
过 计算 水 库 上 游 的 面 雨量 .最 大 汇 水 量 以 及 洪水 峰 
值 , 并 结合 灾后 调查 数据 分 析 对 比 ,将 探讨 水 库 上 
游 山 洪 发 生 的 过 程 以 及 所 产生 汇 水 量 对 水 库 造 成 
的 威胁 ,以 期 加 深 对 山区 小 流域 暴雨 山洪 形成 、 发 
展 和 致 灾 演 变 的 理解 ,为 小 流域 暴雨 山洪 灾害 预警 
和 风险 评 佑 提供 科学 依据 和 技术 参考 。 


1 资料 与 方法 


1.1 研究 区 概况 

射 月 沟 流 域 位 于 新 疆 哈 密 东 北部 山区 ,为 丘陵 
山地 ,植被 覆盖 稀少 ,土壤 蕾 水 能 力 差 ; 射 月 沟 水 库 
上 游 流 域 面积 406 km。 射 月 沟 水 库 距 离 哈 密 市 
115 km ,该 水 库 总 库容 6.70 x 10° ma ,溢洪道 最 大 洪 
峰 流 量 为 380 m.s ,最 大 坝 高 41 m, 大 坝 长 度 
403 m。 

12 数据 资料 

地 面 降水 观测 数据 由 哈密 市 气象 局 提供 ,哈密 
市 伊 州 区 射 月 沟 流 域 及 周边 2018 年 7 月 31 日 01 时 
~14 时 的 14 个 气象 观测 站 的 地 面 观测 逐 时 雨量 数 
据 。 哈 密 市 伊 州 区 防洪 办 和 水 文 局 提供 的 射 月 沟 
水 库 雨 量 站 数据 (2018 年 7 月 31 日 01 时 ~ 10 时 数 
据 ). 射 月 沟 水 库 库容 以 及 水 库 上 游 暴 十 洪水 人 库 、 
Yaa Tee .洪水 漫 坝 省 口 时 间 节 点 等 信息 资料 。 

多 源 融 合 降水 数据 来 自 于 国家 气象 信息 中 心 
(http://www.nmic.en) ,该 数据 资料 是 基于 地 面 自动 
气象 站 与 CMORPH 卫星 反 演 降 水 数据 ,全国 雷 达 拼 
图 定量 估 测 降水 数据 综合 而 成 的 中 国 区 域 融 合 降 
水 分 析 产 品 (CMPAS: CMA Multi-source merged Pre- 
cipitation Analysis System, 简称 :CMPAS) ,时 空 分 辩 
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率 分 别 为 1h、0.05° x 0.05。 ,该 数据 源 的 核心 算法 利 
用 概率 密度 函数 匹配 法 订正 卫星 雷达 降水 的 系统 
偏差 ,采用 贝 叶 斯 模型 平均 法 .最 优 插值 法 以 及 空 
间 降 尺度 技术 形成 三 源 降水 融合 产品 "1。 

基础 地 理 信息 数据 包括 1:50 000 的 数字 高 程 
模型 (digital elevation model, DEM) ,空间 分 辨 率 25 
m, 土 地 利用 数据 采用 2015 年 还 感 解 译 数据 ,空间 
ST HEE 30 mo 
1.3 流域 面 雨量 的 计算 与 获取 

由 图 1 给 出 射 月 沟 流 域 及 其 周边 气象 站 分 布 示 
意图 可 知 , 射 月 沟 流 域内 气象 站 非常 稀 玻 , 仅 有 小 
堡 气象 观测 站 (Y6344) 和 射 月 沟 水 库 雨 量 站 2 个 
站 。 作 为 模拟 洪水 淹没 的 面 雨量 需要 考虑 表征 实 
际 降水 强度 和 分 布 ,为 减 小 面 雨量 的 估算 误差 , 选 
取 了 射 月 沟 流域 及 周边 附近 总 共 15 个 自动 气象 站 
点 降水 资料 ,采用 泰 森 多 边 形 法 (Thiessen) SMEER 
权重 法 (Inverse Distance Weighting, IDW ) 、 兽 通 克 里 
金 法 (Kriging) 等 3 种 空间 插值 法 估算 射 月 沟 流域 面 
雨量 ,同时 通过 多 源 融 合 降水 数据 (CMAPS) 获 取 射 
月 沟 水 库 上 游 范围 内 的 面 雨 量 。 
1.4 模型 参数 设置 

Floodarea 模型 基于 水 动力 方法 ,同时 考虑 了 一 
个 栅 格 周围 8 个 单位 ,对 邻近 栅 格 单元 的 泄 和 人 量 由 
manning-strickler 公式 "计算 ,坡度 由 单元 之 间 的 最 
低 水 位 和 地 形 的 最 高 高 程 之 间 的 差异 所 决定 ,对 每 
一 个 单元 都 进行 计算 ,位 于 对 角 线 的 单元 ,以 不 同 
的 长 度 算 法 来 计算 。 不 同 于 静态 的 洪水 演进 模拟 ， 


表 1 不 同 土地 利用 类 型 对 应 的 Strickler 系数 ” 
Tab.1 Strickler induce of different types 
土地 窗 盖 类 型 ”耕地 草地 水体 ”居民 地 ”荒漠 
Strickler 系 数 18 50 40 25 33 
Manning 系数 0.06 0.02 0.03 0.04 0.03 


位 :m) o 
1.5 累计 汇 水 量 的 计算 

在 特大 暴雨 情况 下 ,相关 研究 表明 地 表 大 量 的 
汇 水 主要 来 源 于 大 气 降水 , 边 坡 的 渗水 影响 较 小 ™， 
此 本 文 的 累计 汇 水 量 仅 考 虑 大 气 降水 。 在 AreGIS 
环境 下 运行 Floodarea 模 型 ,输出 结果 为 表示 浴 没 水 
深 的 栅 格 数据 ,并 且 可 以 得 到 每 个 时 相 的 洪水 淹没 
范围 和 洪水 历时 特征 。 累 计 汇 水 量 的 计算 为 假定 每 
个 单元 格 有 一 个 单位 的 淹没 水 量 ,依据 水 流 方向 依 
次 计算 每 个 单元 格 所 能 积累 的 水 量 , 计 算 公 式 为 : 

VEh rN (1) 

式 中 池 为 汇 水 量 (m ) ;h 为 淹没 深度 (m) 5 anev 为 
DEM 的 空间 分 辨 率 (m)。 汇 水 量 的 大 小 与 洪水 淹没 
深度 直接 相关 ,通过 Floodarea 模 型 得 出 淹没 深度 有 
后 ,可 根据 淹没 深度 与 机 格 单元 大 小 的 乘积 ,计算 
出 每 个 栅 格 的 汇 水 量 , 逐 时 汇 水 量 的 计算 需 相 应 地 
减 去 前 1 小 时 的 汇 水 量 。 
1.6 流域 边界 提取 

流域 地 理 信 息 关键 的 一 步 首 先是 提取 确定 流 
域 的 范围 , 即 射 月 沟 水 库 上 游 的 集 水 面积 。 在 自然 
流域 的 情况 下 ,通过 ArcGIS 扩展 模块 Are Hydro 
Tools 实现 流域 集 水 面积 的 提取 。 整 个 提取 过 程 如 


Floodarea 模型 在 每 个 时 相 的 运行 过 程 , 即 运行 时 间 
与 相应 淹没 范围 和 水 深 均 以 栅 格 形式 呈现 和 储存 ， 
直观 明了 , 便 以 查询 ,人 研究 证 明 用 以 界定 洪水 淹没 
范围 和 预警 洪涝 风险 的 效果 较 好 -””。 由 于 本 文 
选择 研究 区 为 射 月 沟 水 库 上 游 区 ,不 受 堤 坝 、 道 路 、 
建筑 物 等 对 洪水 模拟 过 程 的 影响 ,选用 Floodarea 模 
型 中 的 暴雨 模块 (Rainstorm ) ,该 模块 的 输入 参数 包 
括 逐 时 面 雨 量 数据 Ut SY DEM 粗糙 度 (Strickler 系 
数 ) 和 最 大 水 流 交 换 率 。 该 模块 中 输入 逐 时 面 雨量 
数值 (单位 :mm*h ') ,最 大 交换 率 值 为 1% ,时 间 步 
长 为 1h。 根 据 土 地 利用 类 型 数据 改进 Strickler 系 
数 ,在 ArcGIS 10.2 的 平台 下 利用 Reclassify 工 具 , 对 
地 表 粗 糙 度 重新 赋值 ,得 到 改进 后 的 射 月 沟 流 域 粗 
烙 度 系数 ( 表 1) ,模型 最 终 的 输出 结果 为 从 地 面 开 
始 的 淹没 水 深 , 其 单位 同 高 程 模型 的 单位 相同 ( 单 


F: (1) 原始 DEM 数 据 预 处 理 , 对 流域 洼地 进行 填 
HE. (2) 基于 D8 算 法 生成 流 回 。(3) 累计 流 生成 。 
(4) 流域 出 口 点 生成 。($) 通过 Batch Sub Watershed 
Definition 生成 射 月 沟 水 库 流域 集 水 面积 ,在 生成 的 
Sub Watershed 矢量 文件 属性 表 中 可 得 出 射 月 沟 水 
库 流 域 面积 为 400 km (接近 实测 集 水 面积 406 
km )。 接 着 利用 选取 得 到 的 流域 边界 裁剪 Floodar- 
ea 模型 模拟 所 需 的 DEM 、 土 地 利用 数据 等 。 


2 射 月 沟 流域 面 雨量 与 汇 水 量 估 算 
及 其 洪水 入 库 分 析 
2.1 射 月 沟 水 库 上 游 面 雨量 计算 及 验证 


选取 哈密 市 伊 州 区 射 月 沟 流 域 周边 2018 年 7 
月 31 日 01 时 ~14 时 的 10 个 气象 观测 站 所 采集 的 逐 
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时 雨量 数据 ,分 别 采用 泰 森 多 边 形 法 (Thiessen) JX 
距离 加 权 法 (IDW) 和 普通 克 里 金 法 (Kriging) 估 算 射 
月 沟 水 库 上 游 逐 时 面 雨 量 。 除 了 上 述 空间 插值 所 
得 面 雨量 外 ,将 获取 的 多 源 融 合 降 水 分 析 产 品 
(CMPAS) 也 作为 该 流域 面 雨量 数据 ,4 种 面 雨量 结 
果 见 表 2。 


由 累计 雨量 统计 结果 可 知 :2019 年 7 月 31 日 01 
时 ~14 时 射 月 沟 水 库 上 游 汇 水 区 域 平均 降水 量 , 泰 
森 多 边 形 法 所 得 的 累计 面 雨量 最 大 , 值 为 109.6 
mm。 其 次 是 多 源 融 合 降水 分 析 产 品 ,累计 面 雨 量 为 
100.7 mm ,普通 克 里 金 法 计算 所 得 的 累计 面 雨 量 最 
小 , 仅 为 58.7 mm, 


R2 射 月 沟 水 库 上 游 区 每 小 时 面 雨 量 平均 值 /mm 


Tab.2 Statistical table of area rainfall per hour in the upper reach of Sheyuegou / mm 


项 目 累计 面 雨量 1 时 2 时 3 时 4 时 5 时 6 时 7 时 8 时 9 时 10 时 1 时 12 时 13 时 14 时 
Kriging 58.7 0 1.4 14. 42 34 43 134 128 76 7.3 23 0.1 02 03 
IDW 90.9 0 2.1 1d 8 3.5 4 241 219 13.1 86 3.5 0.2 0.3 0.5 
Thiessen 109.6 0 2.9 15 105 38 42 292 276 161 £88 41 0 04 0.5 
CMPAS 100.7 0 34 24 108 39 54 247 247 143 63 34 06 04 04 


为 了 验证 面 雨量 估算 结果 的 可 靠 性 ,利用 未 参 
与 面 雨 量 估算 的 气象 站 逐 时 降水 数据 (其 中 射 月 沟 
水 库 雨量 站 为 31 日 01 时 ~09 时 ), 提 取 每 站 点 相对 
应 的 逐 时 降雨 估算 值 ,得 出 面 雨量 相关 系数 值 在 
0.76 ~ 0.97 之 间 , 均 通 过 信和 度 0.05 的 检验 ,相关 性 较 


2.2 射 月 沟 水 库 上 游 降 水 汇 水 量 分 析 

为 进一步 讨论 面 雨 量 对 流域 汇 水 量 的 影响 , 基 
F Floodarea 模型 分 别 将 Thiessen IDW , Kriging 和 
CMPAS 的 逐 时 面 雨量 作为 输入 因子 ,通过 累计 汇 水 


好 ( 表 3)。3 种 插值 法 的 估算 结果 与 实况 比较 ,相关 
系数 均 在 0.75 以 上 , 均 方 根 误 差 均 超过 了 2.23; 而 
CMPAS 面 雨量 值 与 实况 值 间 的 相关 系数 以 及 决定 
系数 都 高 于 空间 插值 法 的 佑 算 值 ,表现 最 优 。 


表 3 ” 面 雨 量 估算 值 与 地 面 观 测 值 的 相关 性 分 析 
Tab.3 Correlation analysis between estimations 


and observations 


相关 系数 P 决定 系数 RR ”” 均 方 根 误差 RMSE 
Kriging 0.76 0.58 2.25 
IDW 0.84 0.71 2.23 
Thiessen 0.80 0.64 2.80 
CMPAS 0.97 0.94 1:15 
(a) 
3000) Kuging 
全 Thiessen 


~- CMPAS 


累计 汇 水 量 /10:me 
a 
s 


量 计算 公式 ,得 出 射 月 沟 水 库 上 游 降 水 汇 水 量 与 洪 
峰 流 量变 化 图 (图 3)。 

由 射 月 沟 水 库 上 游 累 计 汇 水 量变 化 图 (图 3 a) 
可 知 : 4 种 面 雨量 计算 所 得 此 次 强 降水 在 射 月 沟 水 
库 上 游 累 计 汇 水 量 为 1.47 x 10 ~ 2.89x10 ms ,达到 
射 月 沟 水 库 总 库容 (6.70 x 10° m) HI 2.2 ~ 4.3 1% A 
洪水 增长 趋势 得 知 ,07 ~ 10 时 为 山洪 暴涨 时 段 ,这 
与 射 月 沟 流 域 自动 南 量 观测 站 的 小 时 暴雨 强度 一 
致 ,3 最 大 汇 水 量 在 1.08x10 ~ 1.86x10 m ÈM- 
从 洪峰 流量 变化 趋势 均 可 以 明显 看 出 (图 3b),7 月 
31 日 02 时 ~14 时 出 现 两 个 峰值 ,第 一 个 山洪 高 峰 
出 现在 05 时 ,第 二 个 峰值 出 现在 08 时 ,其 中 由 
Thiessen 最 大 洪峰 流量 达 2 039 ms, 而 Kriging 最 


(b) 

æ Kriging 
æ IDW 

全 Thiessen 
~- CMPAS 
> 最 大 洪峰 流量 


洪峰 流量 /ma s-! 
js 
© 
© 
© 


图 3 射 月 沟 水 库 上 游 降 水 汇 水 量 与 洪峰 流量 曲线 图 


Fig. 3 Curves of water-collecting amount and the flood peak flow in the upper reaches of Sheyuegou reservoir 
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大 洪峰 流量 为 928 ms …。 根 据 水 利 部 门 灾后 调查 
的 最 大 洪峰 流量 (1 848 m.s- !) 与 CMPAS 最 大 洪峰 
流量 (1 756 m-s ') 相 比 , 两 者 误差 最 小 ,精确 度 达 
到 95%。 通 过 面 雨量 和 最 大 洪峰 量 的 结果 验证 ,多 
源 融 合 降水 产品 (CMPAS) 的 估算 结果 的 可 靠 性 和 
精确 性 表现 最 优 ,由 此 推断 此 次 射 月 沟 水 库 上 游 区 
山洪 累计 汇 水 总 量 达到 2.64x10’ m°. 
2.3 ” 射 月 沟 水 库 上 游 洪 水 入 库 分 析 

图 4 是 基于 CMPAS 的 射 月 沟 水 库 上 游 的 汇 水 
量变 化 过 程 。 射 月 沟 水 库 总 库容 为 6.70 x 10° mm’, Yat 
洪 道 最 大 洪峰 流量 为 380 m-s, 7K ERT LAB 7K 8.07 
x 10° mm, 并 保持 正常 运行 ,7 月 31 日 02 时 水 库 库容 
量 为 4.04 x 10° mm。 通过 累计 汇 水 量 与 水 库 溢洪道 
开始 溢 洪 .洪水 漫 坝 时 间 推 算 , 当 洪 水 入 库 水 量 到 
2.66 x 10 mm 时 , 游 洪 道 开 始 溢 洪 (7 时 35 分 ,调查 时 
间 ,以 下 同 ), 从 累计 汇 水 量 达 到 这 一 量 值 的 时 间 段 
在 7 月 31 日 05 时 左右 ,相差 约 2h 35 min; 当 洪水 人 
库 水 量 到 4.03 x 10° mm 时, 水库 开始 漫 坝 (8 时 左 
右 ) ,而 累计 汇 水 量 达到 这 一 量 值 约 为 06 时 10 分 ， 
相差 1h 50 min。 由 此 推断 ,从 射 月 沟 水 库 上 游 汇 水 
量 形成 山洪 最 后 到 达 水 库 的 时 间 可 能 需要 2 h 左 
右 。 当 水 库 漫 坝 发 生 近 1 h45 min tF, KIR ÆW 
口 (9 时 47 分 ) ,此 时 洪水 超过 大 坝 顶 2 m( 现 场 调查 
数据 ) ,估计 洪水 人 库 量 可 能 达到 了 1.32 x 10’ m° 
值得 注意 是 水 库 大 坝 形成 演 口 后 , 射 月 沟 水 库 上 游 
W5A ZY 1.30 x 107 m HEKA RKE A i UF o 


一 - 累计 汇 水 量 m 逐 时 汇 水 量 
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图 4 基于 CMPAS 的 射 月 沟 水 库 上 游 逐 时 汇 水 量变 化 图 
Fig.4 Based on CMPAS of the curve of the water-collecting 


amount in the upper reach of Sheyuegou reservoir 


3 结论 


x~ 


山地 小 流域 面 雨量 的 估算 ,在 缺少 加 密 降 水 观 


测 站 网 的 情况 下 ,可 利用 小 流域 周边 降水 站 网 数据 
采用 不 同 插值 方法 计算 流域 逐 时 面 雨 量 , 也 可 采用 
国家 气象 信息 中 心 生 产 的 多 源 融 合 降水 逐 时 产 
品 。 通 过 对 4 种 面 雨 量 估算 及 相关 分 析 及 验证 ， 
CMPAS 面 雨量 相关 系数 .决定 系数 为 0.97、0.94 为 
最 大 , 均 方 根 误差 为 1.15 为 最 小 ,CMPAS 面 十 量 优 
于 其 他 3 种 插值 法 ,表现 为 最 优 。 

基于 4 种 面 雨 量 驱 动 Floodarea 模 型 得 出 的 汇 
水 量 在 1.47 x 10’ ~ 2.89 x 10m 之 间 ,最 大 洪峰 流 
量 在 928~2 039 m° -s~ :之 间 , 结 果 差 异 明 显 。 
CMPAS 估算 水 库 上 游 的 最 大 洪峰 流量 1 756 m 
s ,与 水 利 部 门 事后 调查 的 最 大 洪峰 流量 为 1 848 
ms 4H EE ,精确 度 达到 95% ,由 CMPAS 得 到 的 
射 月 沟 水 库 上 游 暴 十 山洪 总 水 量 可 能 为 2.64 x 
10 am? , 远 超 该 水 库 的 防洪 库容 和 溢 洪 承 载 能 

使 用 CMPAS 面 雨 量 计 算 的 汇 水 量 与 灾后 调查 
的 水 库 洪 洪 . 漫 坝 和 溃 口 调查 时 间 进 行 分 析 对 比 ， 
汇 水 量 与 实际 洪水 人 库 有 一 定 的 滞后 性 ,可 能 与 降 
水 强度 分 布 以 及 山顶 汇流 到 达 水 库 的 时 间 等 有 关 ; 
虽然 能 推算 出 射 月 沟 水 库 上 游 的 洪水 致 水 库 淡 口 
的 汇 水 总 量 ,但 由 于 面 十 量 驱 动 floodarea 模 型 的 时 
间 步 长 为 小 时 ,还 不 能 更 为 精准 地 推演 分 钟 时 段 洪 
水 量 影响 水 库 过 程 。 

受 限 于 水 文 资料 的 验证 .以 及 DEM 的 精度 、 产 
流 等 多 方面 因素 影响 ,Floodarea 模 型 模拟 得 到 汇 水 
量 可 能 存在 偏差 ,精准 反 演 小 流域 暴雨 山洪 的 面 十 
量 、 汇 水 量 等 山洪 灾害 防御 关键 技术 是 今后 需要 不 
断 完善 的 工作 。 
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Estimation and application of water-collecting amount during an extreme 
heavy rainfall induced flash flooding in Hami City on 31st July 2018 


WANG Mei-xia, ZHANG Tai-xi, YU Xing-jie, ZHANG Lian-cheng, ZHANG Xu 
( Xinjiang Climate Center , Urumqi 830002 , Xinjiang , China) 


Abstract: In order to mitigate the flooding disasters induced by severe rainfall over mountain regions, it is very 
important to estimate the water collection amount. However, in arid mountain regions , meteorological and hydrologi- 
cal stations are usually scarce and unevenly distributed , leading to difficulties in the estimation of this parameter. In 
such cases, a reasonable estimation of the areal rainfall is necessary and offers the initial conditions for driving hy- 
drological models. In this paper, we studied an extreme rainfall event that occurred in Hami City, Xinjiang, China 
on July 31,2018 and that resulted in the break of a local dam. The northeastern part of Hami corresponds to the 
easternmost area of the Tianshan Mountains. This area is characterized by very hot summers and an annual mean 
precipitation of 43 — 118 mm. Due to a lack of vegetation, flooding disasters occur frequently along the hilly south- 
ern slope of this area; nevertheless, local meteorological and hydrological stations are rare. On July 31,2018, from 
midnight to 13:00, an extreme rainfall event occurred over northeastern Hami and led to a severe flood. This caused 
the break of the Sheyuegou Reservoir and, consequently , the death of many people. In this paper, the areal rainfall, 
water collection amount, and breaking process of the local dam were objectively evaluated. In order to obtain a rea- 
sonable estimation of the areal rainfall, three interpolation methods were applied to the meteorological station data, 
together with a multi-source fusion data reanalysis. Overall, the four estimation methods demonstrated correlations 
between 0.76 — 0.97; notably , the best estimation was given by the multi-source fusion data reanalysis. The estimat- 
ed areal rainfall was then used to drive the FloodArea Model, obtaining the water collection amount for the upper 
reach of the Sheyuegou Reservoir; the values varied between 1.47110’ m° and 2.892x 10’ m°. The water collection 
amount obtained by using the CMPAS areal rainfall was further compared with the results of a post-disaster investi- 
gation on the dam breaking process. The calculated hourly water collection amount and the flood peak flow were in 
good agreement with the results of the post-disaster investigation. In particular, the flood peak flow was 1.756 x 10° 
ms ' (very similar to that derived from hydrological observations ) , while the total flood charge was 2.639x 10’ m° 


(much higher that the flood control and spillway carrying capacities of the reservoir ). 
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